Coquaternionic quantum dynamics for two-level systems 
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The dynamical aspects of a spin-^ particle in Hermitian coquaternionic quantum 
theory is investigated. It is shown that the time evolution exhibits three different 
characteristics, depending on the values of the parameters of the Hamiltonian. When 
energy eigenvalues are real, the evolution is either isomorphic to that of a complex 
Hermitian theory on a spherical state space, or else it remains unitary along an 
open orbit on a hyperbolic state space. When energy eigenvalues form a complex 
conjugate pair, the orbit of the time evolution closes again even though the state 
space is hyperbolic. 
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Over the last decade or so there have been considerable interests in the study of com- 
plexified dynamical systems; both classically [IH6] and quantum mechanically [7HT7]. For 
a classical system, its complex extension typically involves the use of complex phase-space 
variables: (x,p) — > (x + ixi,Po + ipi)- Hence the dimensionality of the phase space, i.e. 
the dynamical degrees of freedom, is doubled, and the Hamiltonian H(x,p) in general also 
becomes complex. For a quantum system, on the other, its complex extension typically 
involves the use of a Hamiltonian that is not Hermitian, whereas the dynamical degrees of 
freedom associated with the space of states — the quantum phase space variables — are kept 
real. However, a fully complexified quantum dynamics, analogous to its classical counter- 
part, can be formulated, where state space variables are also complexified [El 119] . 

The present authors recently observed that there are two natural ways in which quantum 
dynamics can be extended into a fully complex domain [19] , where both the Hamiltonian and 
the state space are complexified. In short, one is to let state space variables and Hamiltonian 
be quaternion valued; the other is to let them coquaternion valued. The former is related 
to quaternionic quantum mechanics of Finkelstein and others [2U1 ETJ, whereas the latter 
possesses spectral structures similar to those of PT-symmetric quantum theory of Bender 
and others [7HT0]. The purpose of this paper is to work out in some detail the dynamics of 
an elementary quantum system of a spin-1 particle under a coquaternionic extension, in a 
manner analogous to the quaternionic case investigated elsewhere [22]. 

As illustrated in [19], a coquaternionic dynamical system arises from the extension of 
the real and the imaginary parts of the state vector in the complex- j direction, where j 
is the second coquaternionic 'imaginary' unit (described below). The general dynamics is 
governed by a coquaternionic Hermitian Hamiltonian, whose eigenvalues are either real or 
else appear as complex conjugate pairs. Here we examine the evolution of the expectation 
values of the five Pauli matrices generated by a generic 2x2 coquaternionic Hermitian 
Hamiltonian. We shall find that, depending on the values of the parameters appearing in 
the Hamiltonian, the dynamics can be classified into three cases: (a) the eigenvalues of 
H are real and the dynamics is strongly unitary in the sense that the 'real part' of the 
dynamics on the reduced state space is indistinguishable from that generated by a standard 
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complex Hermitian Hamiltonian; (b) the eigenvalues of H are real and the states evolve 
unitarily into infinity without forming closed orbits; and (c) the eigenvalues of H form a 
complex-conjugate pair but the dynamics remains weakly unitary in the sense that the real 
part of the dynamics, although generating closed orbits, no longer lies on the state space of a 
standard complex Hermitian system. Interestingly, properties (b) and (c) are in some sense 
interchanged in a typical PT-symmetric Hamiltonian where the orbits of a spin-| system 
are closed when eigenvalues are real and open otherwise. These characteristics are related 
to the three cases investigated recently by Kisil [23] in a more general context of Heisenberg 
algebra, based on the use of: (i) spherical imaginary unit i 2 = — 1; (ii) parabolic imaginary 
unit i 2 = 0; and (iii) hyperbolic imaginary unit i 2 = 1. The use of coquaternionic Hermitian 
Hamiltonians thus provides a concise way of visualising these different aspects of generalised 
quantum theory. 

Before we analyse the dynamics, let us begin by briefly reviewing some properties of 
coquaternions that are relevant to the ensuing discussion. Coquaternions [21] . perhaps 
more commonly known as split quaternions, satisfy the algebraic relation 

i 2 = -1, j 2 = k 2 = ijk = +1 (1) 

and the skew-cyclic relation 

ij = —ji = k, jk = —kj = —i, ki = —ik = j. (2) 

The conjugate of a coquaternion q = q + iq\ + jq 2 + kq 3 is q = qo — iq\ — jq 2 — kq 3 . It 
follows that the squared modulus of a coquaternion is indefinite: qq = q 2 + q 2 — q\ — q\- 
Unlike quaternions, a coquaternion need not have an inverse q~ x = q/{qq) if it is null, i.e. 
if qq = 0. The polar decomposition of a coquaternion is thus more intricate than that of a 
quaternion. If a coquaternion q has the property that qq > and that its imaginary part 
also has a positive norm so that q\ — q\ — q\ > 0, then q can be written in the form 

q= |g|e^ = |g|(cos0 ff + i g sin0 g ), (3) 

where 

iqi+JQ2 + kq 3 x 
i q = — and U q = tan 

That a coquaternion with 'time-like' imaginary part admits the representation ^ leads to 
the strong unitary dynamics generated by a coquaternionic Hermitian Hamiltonian. On the 
other hand, if qq > but q\ — q\ — q\ < 0, i.e. if the imaginary part of q is 'space-like', then 

q= \q\e iqeq = \q\(cosh9 g + i q smhO g ), (5) 

where 

i q = i ^+M±iM and 6 q = tanh- 1 ( ^ + f + ?i 

V-ii + 4 + 4 V |go1 

If qq > and 4 ~ 4 ~ 4 = 0> then q = q (l + i q ), where i q = q$ (iqi + jq 2 + kq 3 ) is the 
null pure-imaginary coquaternion. Finally, if qq < 0, then we have 

q= |g|e i9e9 = \q\(smh 9 q + i q cosh 9 q ), (7) 
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where 

+ i + i V * / 

As indicated above, the fact that the polar decomposition of a coquaternion is represented 
either in terms of trigonometric functions or in terms of hyperbolic functions manifest itself 
in the intricate mixture of spherical and hyperbolic geometries associated with the state 
space of a spin-| system, as we shall describe in what follows. 

In the case of a coquaternionic matrix H, its Hermitian conjugate H^ is defined in a man- 
ner identical to a complex matrix, i.e. W is the coquaternionic conjugate of the transpose 
of H. Therefore, for a coquaternionic two-level system, a generic Hermitian Hamiltonian 
satisfying W = H can be expressed in the form 

5 

H = u t + ^uia u (9) 



i=i 



1 \ . / -i \ A / 1 



where {wi}z=o..5 £ H, and 

*i = 1 1 o), ** = o ;» a *-yo -i 

are the coquaternionic Pauli matrices. The eigenvalues of the Hamiltonian ^ are given by 



E± = U Q ± y u l + ^2 + U \ ~ U l ~~ U 5- (ii) 

Thus, they are both real if u\ + u\ + u\ > ul + u\; otherwise they form a complex conjugate 
pair. This, of course, is a characteristic feature of a PT-symmetric Hamiltonian. 

A unitary time evolution in a coquaternionic quantum theory is generated by a one- 
parameter family of unitary operators e~ At , where A is skew-Hermitian: = —A. As 
in the case of complex quantum theory, we would like to let the Hamiltonian H be the 
generator of the dynamics. For this purpose, let us write 

i — - (iu 2 + ju 4 + ku 5 ), (12) 



where v = \/u\ — u\ — u\ if u\ + u\ < ul, and v = \/ u\ + u\ — u\ if u\ < u\ + u\. Then we 
set A = iH and the Schrodinger equation in units h = 1 is thus given by (cf. [22] ) 

\f) = -%H\f). (13) 

It is worth remarking that when u\-\-u\ < u\ we have i 2 = —1, whereas when u\ < u\-\-u\ we 
have i 2 = +1. In either case iH is a skew-Hermitian operator satisfying {iH)^ = —iH; thus 
Q-iHt f orma Uy generates a unitary time evolution that preserves the norm = ipiipi + 

where is the coquaternionic conjugate of ip so that represents the Hermitian 
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conjugate of The conservation of the norm can be checked directly by use of the explicit 
form of the Schrodinger equation in terms of the components of the state vector: 



i>i\ = ( ~{uo + u 3 )iif)i - uiiip 2 - v^i 
4>2 J \-(uo-u 3 )i^2-uii^i + vifi 



(14) 



Here we have assumed u\ + u\ < u\ so that v — \fu\ — u\ — u\] if w| < u\ + u\, we have 
v = \Ju\ + ul — u\ and the sign of v in ( 14 ) changes. 



To investigate properties of the unitary dynamics generated by the Hamiltonian (|9j) we 
shall derive the evolution equation satisfied by what one might call a 'coquaternionic Bloch 
vector' a, whose components are given by 

(ip\ai\ip) 

°i = -(^p / = !,•• -,5. (15) 



By differentiating (j\ in t for each I and using the dynamical equation (14), we deduce, after 



rearrangements of terms, the following set of generalised Bloch equations: 



1 • ^3 / . x 

Hffl = VOj, \U 2 (J 2 + M 4 0"4 + U^Os) 

1 V 

i • !/ \ 
^a 2 = -{u 2 u 3 ai - uiu 2 a 3 + u u^a A - uqu^) 

z v 

\a 3 = —v&\ H — -{u 2 a 2 + u 4 a 4 + u 5 cr 5 ) (16) 

i- !/ \ 
~o"4 = -(— u 3 U4(Ji + u u 5 a 2 + UiU A a 3 + Mono's) 

i • !/ x 

where we have assumed u\ + u| < u\ so that v = \/u^ — u\ — u\. This is the region in 
the parameter space where the coquaternion appearing in the Hamiltonian has a time-like 
imaginary part. Note that these evolution equations preserve the condition: 

a\ + a\ + a\ - a\ - a\ = 1, (17) 

which can be viewed as the defining equation for the hyperbolic state space of a coquater- 
nionic two-level system. 

Let us now show how the dynamics can be reduced to three-dimensions so as to provide 
a more intuitive understanding. For this purpose, we define the three reduced spin variables 

<J x = (Tx, cr y =^(u 2 a 2 + u i a 4: + u 5 a 5 ), a z = a 3 . (18) 
We can think of the space spanned by these reduced spin variables as representing the 'real 



part' of the state space (17). Then a short calculation making use of (16) shows that 

i 



2 a x = vcr z - u 3 a y 
\(T y = u 3 o x - u x o z (19) 
\o z = u x o y - ua x , 



FIG. 1: (colour online) Dynamical trajectories on the reduced state spaces. In the parameter region 
u\ > u\ + u\ the reduced state space is just a two-sphere, upon which the dynamical equations 
(19) generate Rabi oscillations (left figure). In the parameter region u\ < u\ + it§ the reduced 
state space is a two-dimensional hyperboloid, and the dynamical equations (26) generate open 
trajectories on this hyperbolic state space, if the energy eigenvalues are real (right figure). If the 
eigenvalues are complex, the open trajectories are rotated to form hyperbolic Rabi oscillations. 



or, more concisely, a - 
coquaternionic spin-^ 



= 2B x a where B = (u\, v, U3 
system is a hyperboloid 



Hence although the state space of a 



17|), remarkably in the region u\ + u h 



< ui 



the reduced spin variables a x , a y , a z denned by ( |18p obey the standard Bloch equations (19). 
In particular, the reduced motions are confined to the two sphere S 2 : 



a l + a v + a l = const., 



(20) 



where the value of the right side of (20) depends on the initial condition (but is positive 



and is time independent). Put the matter differently, in the parameter region u A + ui, < 
u 2) the dynamics on the reduced state space S 2 induced by a coquaternionic Hermitian 
Hamiltonian is indistinguishable from the conventional unitary dynamics generated by a 
complex Hermitian Hamiltonian. This corresponds to the situation in a PT-symmetric 
quantum theory whereby in some regions of the parameter space the Hamiltonian is complex 
Hermitian (e.g., a harmonic oscillator in the Bender-Boettcher Hamiltonian family H = 
p 2 +x 2 (ix) e [7] , or the six-parameter 2x2 matrix family in [25J ) . Some examples of dynamical 
trajectories are sketched in figure [TJ 

The evolution of the other dynamical variables a 2 , 04, cr 5 can be analysed as follows. Recall 



that the dynamics (19) preserves the relation (20). Thus, by subtracting (20) from (17) and 



rearranging terms we deduce that 



- (M2CT4 + W4O2) 2 + ( M 40"5 - W5O4) 2 - (u 5 a 2 + u 2 a 5 ) 2 = const. 



(21) 



This shows that the evolution of the vector (a 2 , 04, cx 5 ) is confined to a hyperbolic cylinder. 
It turns out that the time evolution of these 'hidden' dynamical variables a 2 , 04, 05 can also 
be represented in a form similar to Bloch equations if we transform the variables according 
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to a yi = W4O5 — 115(74, o V2 = W5O2 + and a y3 = W2O4 + W4O2. In terms of these auxiliary 
variables we have 

1 • u o , , x 

2^1 = - — \U^,(Jy 2 +U i ay i ) 

- — {U2(Ty 3 +U 5 (7y 1 ) (22) 



1 • 

2 vz 



1 . «0 / 



It should be evident that these dynamics are confined to a hyperboloid: 

- a l + °m + ^3 = const - ( 23 ) 
Note, however, that when u = we have & yi = b y2 = d„ 3 = from (22), while o 2 ,o 4 ,o 5 



are in general evolving in time. Hence in transforming the variables into a yi ,a y2 ,o yz , part 
of the information concerning the dynamics is lost. 



We see from (20) and (21) that on the 'imaginary part' of the state space the dynamics 
is endowed with hyperbolic characteristics, which nevertheless is not visible on the reduced 
state space, or the 'real part' of the state space S 2 spanned by a x , a y , a z . 

When u\ = u\ + u\ so that the imaginary part of the coquaternion appearing in the 
Hamiltonian is null, a calculation shows that the reduced spin variables obey the following 
dynamical equations: 



\a x = -u 3 a y 

1 j. 

2 

\o z = UxCTy, 



v y = ~u 3 a x + u\a z (24) 



and preserve o 2 x — a 2 + a 2 z . 

When u\ < u\ + u\ so that the imaginary part of the coquaternion in the Hamiltonian 
is space-like, the structure of the state space, as well as the dynamics, change, and they 
exhibit an interesting and nontrivial behaviour. The five-dimensional spin variables in this 
case evolve according to 



2 



1 u 3 \ 

iCTi = -ZAX3 (W 2 0"2 + W4O4 + M5O5) 

V 

1 • !, \ 

2O2 = -(W2W3C1 - «1«20"3 + u u 5 a 4 - U U 4 (T 5 ) 

|<7 3 = Z/Oi + — (u 2 2 + M40" 4 + «50"5) (25) 

l- !/ \ 

^0 4 = -(-U3M4O1 + U U 5 a 2 + M1M4O3 + U Q U 2 (T 5 ) 
z V 

1- l , X 

o^5 = -(-M3M5O1 - UqU^2 + UiU 5 a 3 - U U 2 <J4), 



where v = \fu\-\-u\ — u 2 ,. These evolution equations preserve the normalisation (17). 
However, in the region u 2 < u\ + u\ the reduced spin variables a x , a y , a z defined by (18) no 



longer obey the standard Bloch equations (19); instead, they satisfy 

i • 



2-x — ~vo z - u 3 a y 
l&y = -u 3 cr x + u x a z (26) 
\a z = u x o y + va x , 
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FIG. 2: (colour online) Conic sections and PT phase transition: changes of orbit structures. A 
projection of the orbits on the hyperboloid, for parameters just above the transition to complex 
energy eigenvalues, is shown on the left side. The orbits form circular sections. On the right side 
we plot orbits of hyperbolic Rabi oscillations further into complex energy eigenvalues. The energy 
eigenvalues determine the angle between the axis of rotation and the axis of the hyperboloid. When 
eigenvalues are complex, the axis of rotation is within the hyperboloid, leading to closed orbits on 
the state space generated by circular sections. When the imaginary part of the coquaternion 
appearing in the Hamiltonian is null, we have parabolic sections of the hyperboloid; whereas when 
the energy eigenvalues are real, the angle of the two axes is larger than 7r/4, and open orbits are 
generated by hyperbolic sections. 



and preserve the relation 

ol - °l + °l = const. (27) 
We thus see that at the level of reduced spin variables in three dimensions, the state space 



changes from a two-sphere (20) to a hyperboloid (27), as the parameters u 2 , «4, u 5 appearing 
in the Hamiltonian change. This transition corresponds to the transition from a complex 
Hermitian Hamiltonian into a PT-symmetric non-Hermitian Hamiltonian. Since the energy 
eigenvalues can still be real even when u\ <u\ + u\, we expect the dynamics to exhibit two 
distinct characteristics depending on whether the reality condition u\ + 11% + u\ > u\ + u\ is 
satisfied. Indeed, we found that on a hyperbolic state space, orbits of the unitary dynamics 
associated with real energies are the ones that are open and run off to infinities. Conversely, 
when the reality condition is violated, these open orbits are in effect Wick rotated to generate 
closed orbits. These features can be identified by a closer inspection on the structure of the 
underlying state space, upon which the dynamical orbits lie. In particular, (|26l shows that 



the dynamics generates a rotation around the axis (u\,u, U3); whereas the state space (27) 



is a hyperboloid about the axis (0,1,0). We have sketched in figure dynamical orbits 



resulting from (26), indicating that there indeed is a transition from open to closed orbits 
as real eigenvalues turn into complex conjugate pairs. 

Intuitively, one might have expected an opposite transition since in a PT-symmetric model 
of a spin- \ system the renormalised Bloch vectors on a spherical state space follow closed 
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orbits when eigenvalues are real, whereas sinks and sources are created when eigenvalues 
become complex [15]. The apparent opposite behaviour seen here is presumably to do 
with the facts that the underlying state space is hyperbolic, not spherical, and that no 
renormalisation is performed here. In figure [2] we have sketched some dynamical trajectories 
when energy eigenvalues are complex. A projection of the dynamical orbits from the a z axis 
(for the choice of parameters used in these plots) shows in which way the topology of the 
orbits are affected by the reality of the energy eigenvalues. 

The evolutions of the other dynamical variables a 2 , 04, 05 are confined to the space char- 
acterised by the relation 

(u 2 <Ja + u 4 a 2 ) 2 - - M5C4) 2 + (u 5 a 2 + u 2 a 5 ) 2 = const., (28) 



instead of the relation (21) of the previous case. However, if we define, as before, three 



auxiliary variables a yi = u^a^ — u^a^ a V2 = u^a 2 + u 2 a^, and a V3 = u 2 o± + u^a 2 , then the 



dynamical equations satisfied by these variables are identical to those in (22), except, of 
course, that the definition of v is different. 

It is interesting to remark that when the imaginary part of the coquaternion appearing in 
the Hamiltonian is space-like, the imaginary unit i has the characteristic of a 'double number' 
or a 'Study number' introduced by Clifford [25], that is, i 2 — 1. Quantum theories generated 
by such a number field (instead of the field of complex numbers) and other hyperbolic 
generalisations, as well as various issues that might arise from such generalisations, have 
been discussed by various authors (e.g., [271 EH]; see also [23] and references cited therein). 
The use of coquaternionic Hermitian Hamiltonian thus captures dynamical behaviours of 
different generalisations of quantum mechanics in a simple unified scheme. 

We thank the participants of the international conference on Analytic and Algebraic 
Methods in Physics VII, Prague, 2011, for stimulating discussions. EMG is supported by 
an Imperial College Junior Research Fellowship. 
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